Multipair DC-Josephson Resonances in a biased all-superconducting Bijunction 
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An all-superconducting bijunction consists of a central superconductor contacted to two lateral 
superconductors, such that non-local crossed Andreev reflection is operating. Then new correlated 
transport channels for the Cooper pairs appear in addition to those of separated conventional Joseph- 
son junctions. We study this system in a configuration where the superconductors are connected 
through gate-controllable quantum dots. Multipair phase-coherent resonances and phase-dependent 
multiple Andreev reflections are both obtained when the voltages of the lateral superconductors are 
commensurate, and they add to the usual local dissipative transport due to quasiparticles. The 
two-pair resonance (quartets) as well as some other higher order multipair resonances are 7r-shifted 
at low voltage. Dot control can be used to dramatically enhance the multipair current when the 
voltages are resonant with the dot levels. 

PACS numbers: 74.50. -|-r, 74.45.-l-c, 73.21. La, 74.78.Na 



One of the most striking manifestation of macroscopic 
quantum coherence is the Josephson effect a DC 
current flows when a phase difference is imposed on a 
junction bridging two superconductors with a narrow in- 
sulating, metallic or semiconducting region. When ap- 
plying a constant voltage bias to this same junction, an 
oscillatory current arises Q and the application of an rf- 
irradiation leads to the observation of Shapiro steps with 
zero differential resistance 0, Q and phase coherence Q ■ 
More generally, the microscopic origin of these effects is 
Andreev reflections of electrons and holes at the bound- 
aries of the two superconductors. The same mechanism 
participates in the appearance of a subgap structure in 
highly transparent voltage-biased junctions, a feature un- 
derstood to be due to dissipative quasiparticle emissions 
called multiple Andreev reflections (MAR) (^,0], which 
were observed in atomic point contact experiments 

Non-local quantum mechanical phenomena Q and en- 
tanglement are nowadays investigated in condensed mat- 
ter physics, in particular in superconducting circuits 0. 
Multiterminal superconducting hybrid devices with one 
superconducting arm and two normal metal electrodes 
have also been studied in the last decade [13] with the 
aim of detecting non-local entangled electron pairs [Tlj . 
There is now convincing experimental data on non-local 
current and noise detection which points in this direc- 
Yet, there is also a growing interest in three- 



12|. 



tion ^ 

terminal all-superconducting hybrid structures [l3l - [l6l |. 
so far mainly in regimes dominated by phase-insensitive 
processes. A recent calculation for a SNS junction, where 
the N region is tunnel-coupled to another superconduc- 
tor, also showed resonances ascribed to voltage-induced 
Shapiro steps [l3|- 



The present work shows that a non-dissipative phase- 
coherent Josephson signal of Cooper pair transport could 
be observed in a device consisting of three superconduc- 
tors driven out of equilibrium. This effect relies on a 
combination of both direct Andreev reflections and non- 
local crossed Andreev reflections (CAR) [13] , and is thus 
directly tied to non-local entangled electron processes 
as well as Josephson physics. Here, the "bijunction" 
which we propose consists of a central superconductor 
S'o coupled via two adjustable quantum dots to two lat- 
eral superconductors Sa and Sh, biased at voltages T4 
and Vb (Fig. [ij. As the coherence length of S'o (which 
is grounded at Vb = 0) is assumed to be larger than its 
width between the dots, this bijunction cannot be simply 
considered as two separated junctions in parallel. Each 
junction consists of a quantum dot, made with e.g. car- 
bon nanotubes l^, l3| or nanowires [13] , and labeled Da 
{a ~ a,b). The dots introduce additional degrees of free- 
dom (position of energy levels, coupling widths) which 
provide full control of the junctions. Equilibrium calcu- 



lations j21| in a similar three-terminal device involving 
normal-metal interfaces showed that a bijunction could 
be a source of spatially correlated pairs of Cooper pairs 
(referred to as "non-local quartets") transmitted into Sa 
and Sh simultaneously. 

This Letter reports on calculations of out-of- 
equilibrium transport in a biased SaDaSoDhSb bijunc- 
tion, with as main results: (i) At commensurate volt- 
ages nVa + mVb = (m and n integers), the two AC 
Josephson terms associated with each bias synchronize 
as DC Josephson resonances, which correspond to the 
phase-coherent transport of n pairs to Sa, and m pairs 
to Sb, from S'o. (ii) The Josephson current-phase rela- 
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FIG. 1: (color online) A Josephson bijunction (left). Super- 
conductors Sa {ct — a, b) are biased at voltages Va, while 5*0 is 
grounded. The distance between the two quantum dot junc- 
tions is comparable to the coherence length. The right panel 
shows the energy diagram for the SaDaSoDi,Si, bijunction and 
a higher order diagram associated with a "sextet" current 
with 3 pairs emitted, 2 in Sa and 1 in 5*6, with Va = —14/2. 



tion of quartet resonances (n = m = 1) and that of some 
higher-order resonances are 7r-shifted at low bias. This 
new mechanism for producing a 7r-shift is of particular 
importance for future interferometry experiments, (iii) 
Gate and/or bias voltages can be tuned to enhance the 
multipair resonances by orders of magnitude as compared 
to the adiabatic regime, making them easily observable 
in experiments, (iv) A DC quasiparticle-pair interference 
term, corresponding to phase-dependent MAR, emerges 
from the dissipative Josephson component. 



A simple phase argument (2l| suggests the existence 
of quartet resonances in a bijunction. Starting with an 
equilibrium situation, the current-phase relation of a sin- 
gle tunnel junction SaSo with phases ipa in Sa and cpo 
in 5*0 is Ic sin {ipa — <Po), to which higher-order harmon- 
ics can also contribute. In a SaSoSb bijunction (with 
phase (pb in Sb)-, there exists in addition a quartet and 
a pair cotunneling supercurrent. The DC quartet su- 
percurrent can be viewed as a non-local second-order 
harmonic Iq = /go sin((pa + Vb ^ 2(po), while the pair 
cotunneling corresponds to a DC Josephson effect be- 
tween Sa and Sb through Sq: Ipc = IpcosiD.{ipa — fb) 
2l|. More generally, assuming large enough trans- 
parencies, multipair currents la/b in electrodes Sa/Sb 
are obtained when differentiating the Josephson free en- 
ergy with respect to the superconducting phases (as- 
suming (/?o = 0) la/b = T,n,m ^a/b,in,m) Sm{n(pa + nUfb) . 

When voltages Va/b are applied to Sa/b, fa and cpb ac- 
quire a time dependence and in the special case where 
nVa + mVb = 0, the adiabatic approximation yields 
d{ntpa{t) + mipb(t)) /dt ~ 0. The corresponding current 
component Ia/b,{n,m) sin{'rnpa{t) + rrnpb(t)) and its higher 
harmonics are constant in time despite the applied volt- 
ages, thus leading to a DC current signaling the existence 
of a multipair resonance. An example of such a reso- 
nance is provided in Fig. [T] from a diagrammatic point 
of view, showing the case 214 + Vb = to lowest or- 
der. The voltage constraint allows to close a resonance 
path provided by one Andreev reflection in Sb and 5*0, 



two in Sa as well as two CAR amplitudes in Sq. Note 
that in general these multipair resonances must coexist 
with the usual AC components /o,(i.o) = loa sin(pa{t) and 
h,{o,i) = Job sin ipbjt), as well as with the MAR currents 
discussed in Ref. [IJ] . While this low-bias argument sug- 
gests the possibility of multipair resonances also at higher 
voltages, an exact out-of-equilibrium calculation at arbi- 
trary voltage is still lacking, and it is discussed below. 

The model Hamiltonian of the SaDaSgDbSb bijunction 
is written a.s H ~ J^j iij + Hd + T^t where T-Lj is the 
Hamiltonian for the lead Sj (j = 0, a, &), expressed with 
the Nambu spinors 



„ (1) 

with Pauli matrices acting in Nambu space. "Hd = 
Y.s,a=a,b^adl^sd'as describcs the single energy level of 
each dot, while the tunneling between the dots and the 
electrodes reads 



nrit) 



jka 
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where dj^ = (dj^^ dai) is the Nambu spinor for electrons 
on dot a. The phases are specified by the applied volt- 
ages ^Pj{t) = (pj"^ + 2eVjt/h. The superconducting gaps 
A are assumed identical and the couplings are assumed to 
be symmetric. The "bare" phases ff\ which are usually 
unimportant in an out-of-equilibrium setup, are relevant 
here in the transport calculations. The width of the su- 
perconducting region Sq is assumed to be negligible, a 
situation which corresponds to the maximum coupling 
between the two junctions forming the bijunction. 

We use a Keldysh path-integral technique, where the 
quadratic electronic degrees of freedom from all three 
leads can be integrated out. The Green's function G 
of the dots which is norrjjerturbative in Tir, is obtained 
from a Dyson equation [Gl involving the free dots Green's 
function, and electrode self-energies with both local and 
non-local propagators (see 23|, and the EPAPS for de- 
tails). The average current from electrode j can then 
be computed using a Meir-Wingreen type formula 
generalized to superconductors [23|, |2^ : 



(/,(t)) = -2Re tr 



{t,t) 



(3) 



where tr stands for the trace in the Nambu-dots space, 
is the self-energy of electrode j, refer to the Keldysh 
component, and o stands for a convolution in time. For 
arbitrary voltages Va and Vb, the time-dependence of the 
system is described in terms of two independent Joseph- 
son frequencies uja = 2eVa/h and Wb = 2eVb/h, the 
Green's function G{t, t') is a function of two times, and 
solving the Dyson equation is a daunting task. However, 
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FIG. 2: (color online) "Broad" dots regime (metallic junc- 
tions): |£a,6| = 6A, Ta,b = 4A. Left: Amplitudes of the 
phase-dependent DC current {lai'pl^'')) main reso- 

nances (with |n| + |m| < 3), in units ofeA/h, centered around 
the values of the phase-independent current (small horizon- 
tal bars). Horizontal axis is Va/Vb, with Vf,/A = 0.3. Upper 
right: Current (la), for the resonance 2Va + VJ, = 0, as a func- 
tion of the phases ip'^^ and ipl^\ showing the dependence in 
2ip'^^ +(p^f^^\ Lower right: the current-phase relation {7a(v'6°')) 
at ifi'^^ = for the resonance Va + Vt 
TT-phase behavior. 



0, which shows the 



when the voltages are commensurate, nVa + mVb = 
(n, m integers), the system becomes time-periodic, and 
the Dyson equation can be solved in Fourier space with 
inversion of large matrices. 

We first consider the regime in which each dot mimics a 
metallic junction, achieved by placing energy levels out of 
resonance > A and choosing large couplings Fq > A 
{a = a,b). Wc compute the DC currents {la/t) for dif- 
ferent ratio of the voltages, satisfying nVa + mVb = 0. 
The results for the largest resonances {\ti\ + \m\ < 3 in 
nVa+mVb = 0) are shown in the left panel of Fig. [2] One 
clearly sees that the resonances are easily distinguished 
from the phase-independent background current. The 
resonant multipair DC-current (/q^^) is a function of the 
combination m^i"^ -I- mipf'\ which implies a simultane- 
ous crossing of n pairs from 5*0 to Sa and m pairs from 
So to Sb- The upper right panel in Fig. [2] shows as an 
example the phase dependence for n = 2 and m ~ 1, 
which is indeed a sinusoidal function of the combination 
Stysi^"* + Vb''^- The existence of DC phase-coherent reso- 
nances despite large nonzero voltages is the result of an 
intrinsic synchronization of the AC currents from each 
junction, which works well beyond the adiabatic regime. 

One of the lowest-order (and larger) resonances cor- 
responds to quartets (T4 = — Vf,), i.e. to the correlated 
transmission of two pairs from 5*0 to Sa and Sb respec- 
tively. The "dual" lowest-order resonance corresponds to 
Vb = Va, where pairs cross from Sa to Sb by cotunncl- 
ing through S'o- The sign of the multipair resonances is 
non-trivial. In particular, the quartet resonance is neg- 



ative, which means that the current {lai^a'^ + Vb'')) is 



of TT-type, as shown in the lower right panel of FiglS] 
Similar sign changes of the multipair current-phase rela- 
tion are also obtained for certain high-order resonances. 
The TT-shift is understood from a simple argument. It 
is related to the internal structure of a Cooper pair via 
the antisymmetry of its wavefunction, similarly to the 
TT-junction behavior of a magnetic junction formed by a 
quantum dot with a localized spin 25|. Starting from 
two Cooper pairs in Sq, the production of a non-local 
quartet consists in forming two non-locally entangled sin- 
glets in the dots Da and Db- These two split pairs cor- 
respond to two CAR amplitudes, as those apparent in 
Figure [H A non-local singlet is obtained by the operator 
71 (4t4; - 4;4t) acting on the empty dots. Applying 
this operator twice to describe a non-local quartet state 
leads to ^Q,D„,L>b = — | ti)a | i't)b, which is recast as the 
opposite of the product of a pair in Da and another one 
in Db- A similar reasoning can be applied in order to 
explain the anomalous sign of higher-order harmonics. 

We now investigate the possibility for optimizing the 
multipair resonances by tuning the dot levels, with Sa = 
—Eb = — 0.4A inside the gap, choosing small values of 
the couplings F^ = Fh = 0.1 A. We focus on the quar- 
tet resonance Va = —Vb for specificity (similar behavior 
is observed for the other resonances). When the bias is 
small enough (Vb ^ 0.1 A here), the system is in the adi- 
abatic regime, and the current does not change when Vb 
is varied. We independently checked with a Matsubara 
formalism calculation (not shown) that this current is the 
same as the one obtained here at equilibrium (Vb = 0). 

The current-phase relations {laif'b'')) ^nd (/f,((y9^°'')) for 
V = 0.09A are shown in the first panel of Fig. [3l These 
average currents are identical, thus are made only from a 
quartet component. They show a purely harmonic func- 
tion of the phase ip^^"^ , and suggest a 7r-junction behavior 
for the quartet resonance near equilibrium. When Vb in- 
creases and the non-adiabatic regime is reached, drastic 
changes appear in the current-phase relations, as shown 
in the next panels of Fig. [31 There, both the sign and 
the (non-sinusoidal) shape of the current-phase relation 
changes rapidly with Vb as it approaches the dot energy 
|£b|. The amplitude of the quartet current near the res- 
onance is ^ 1000 times larger than the one in the adi- 
abatic regime. This resonant effect of the dot levels is 
most apparent by plotting the critical current (the 
maximum of the absolute value of the phase-dependent 
part of lai^p^'^)) as a function of Vb- This is shown in 
the left panel of Fig. [H sharply increases and reaches 
a maximum around Vb — Sb- The large increase in the 
quartet current is due to a double resonant effect: first, 
as the dots have opposite energies Eq = —Sb, the forma- 
tion of a quartet in the double dot as a pair in Da and 
a pair in Db is resonant (this is true for any voltage Vb); 
second when Vb — Sb the tunneling of a pair from Da to 
Sa, and from Db to Sb, is also resonant. 
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FIG. 3: (color online) Current-phase relations /a(</'[°^) (red, 
full curve) and Ii,{tp'^^) (blue, dashed curve), in units of eA/h, 
in the quartet configuration Va = — VJ,, for the resonant dots 
regime: Sa ~ —£b ~ 0.4A and F = 0.1 A. Note the different y 
scales in the different panels. 



Increasing Vf, further, e.g. V& ^ 2 A/3 in the present 
case, we see from the lower panels of Fig. [3] that the 
currents (lai^b^^)) ^^'^ (hiy^'f^)) start to deviate sub- 
stantially. This implies the existence of another phase- 
sensitive process different from the one responsible for 
multipair resonances. We call this current contribution 
jphMAR ^Jqj. phase-sensitive MAR), as it is the result of 
the combination of a multipair process with MAR. The 
lowest order diagram contributing to Ip'^^^R is shown 
in the right panel of Fig. H) It can be seen as the inter- 
ference of the amplitudes of two MAR processes at the 
505*0 and 5(,5o interfaces, each promoting a quasiparticle 
from an energy ~ — A in superconductor Sh to an energy 
^ +A in superconductor Sq. This diagram has a thresh- 
old at V = 2A/3, corresponding to the observed value at 
which /p''^^^^ becomes noticeable. However, unlike the 
usual MAR processes found in single junctions, this pro- 
cess (and similar ones of higher order) has the striking 
property of being phase-dependent. 

Addressing more general resonances, the total DC cur- 
rent can be decomposed into 3 components, as inspired 
by Josephson's work Defining I = {{la), (h)), ffi ~ 
{(fia,'Pb),_V = {Va,Vb), e ={ea,eb) one has J(<^, 1/, e) = 
7^^((^, V, e) + IphMAR^^^ _^ ^qp^y^ Assuming 

electron-hole symmetry to hold, for instance with flat 
normal metal density of states in the leads, one can show 
that in the situation studied here, the DC-current obeys 
the relation: 
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FIG. 4: (color online) Left panel: critical current in units 
of eA/?t as a function of the voltage in the quartet config- 
uration, for resonant dots (same parameters as in Fig. [3|. 
Right panel: Lowest order diagram contributing to the phase- 
dependent MAR current. 



Then the following properties hold: (i) The pure quasi- 
particle current is phase-insensitive and odd in volt- 
ages, (ii) The coherent multipair current I^^^ is a func- 
tion of rup^a^ + mtp^^\ it is odd in phases and even in 
voltages, just like the non-dissipative Josephson term. It 
satisfies m{Ia) = n{Ib)- (hi) The component jp^mar jg 
even in phases and odd in voltages, like the dissipative 
("cos(/3") Josephson component, hut it becomes DC in a 
bijunction. This Jp^^^^^ component is also a function of 
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From an experimental standpoint, multipair reso- 
nances can be directly detected by transport measure- 
ments where one probes the nonlocal conductance ^j^^ 
as a function of Vq, Vf,. The phase coherence of the multi- 
pair current, and its actual dependence in V'i/bi however, 
are more difficult to probe directly. One way would be 
to design specific SQUID geometries or microwave reflec- 
tivity experiments. [27| 

We have shown by non-perturbative out-of-equilibrium 
calculations that coherent multipair and phase- 
dependent MAR processes appear in a superconducting 
bijunction. These are due to crossed Andreev reflection 
processes, through the formation of several entangled 
non-local pairs, and lead to signatures in the DC current 
with very speciflc phase and voltage dependence. A 
natural extension of the present work should focus on 
the role of local Coulomb interaction on the dots [26j . 
While in the resonant case, a simple self-consistent 
mean-fleld treatment should give very similar results, 
away from resonance, the present calculation only works 
qualitatively and a more careful examination is required 
to account for the Coulomb blockade and Kondo regimes. 
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= 



(0) 



(4) 



[1] M. Tinkham, Introduction to Superconductivity (Mc 
Graw-Hill, 1996, Singapore). 



5 



[2] B.D. Josephson, Physics Letters 1, 251 (1962). 

[3] S. Shapiro, Phys. Rev. Lett. 11, 80 (1963). 

[4] C. Vanneste, C. C. Chi, W. J. Gallagher, A. W. Klein- 

sasser, S. I. Raider, and R. L. Sandstrom, J. Appl. Phys. 

64 242 (1988). 

[5] M. Octavio, M. Tinkham, G. E. Blonder, and T. M. Klap- 

wijk, Phys. Rev. B 27 6739 (1983). 
[6] J. C. Cuevas, A. Marti'n-Rodero, and A. Levy-Yeyati, 

Phys. Rev. B54, 7366 (1996). 
[7] E. Scheer, P. Joyez, D. Esteve, C. Urbina, and M. De- 

voret, Phys. Rev. Lett. 78, 3535 (1997). 
[8] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. Lett. 

47, 777 (1935). 
[9] M. Ansmann, H. Wang, R. C. Bialczak, M. Hofheinz, E. 

Lucero, M. Neeley, A. D. O'Connell, D. Sank, M. Weides, 

J. Wenner, A. N. Cleland, and J. M. Martinis, Nature 

461, 504 (2009); L. DiCarlo, M. D. Reed, L. Sun, B. R. 

Johnson, J. M. Chow, J. M. Gambetta, L. Frunzio, S. M. 

Girvin, M. H. Devoret, and R. J. Schoelkopf Nature 467, 

574 (2010). L. DiCarlo, J. M. Chow, J. M. Gambetta, 

Lev S. Bishop, B. R. Johnson, D. I. Schuster, J. Majer, 

A. Blais, L. Frunzio, S. M. Girvin, and R. J. Schoelkopf. 

Nature 460, 240 (2009); 
[10] J. M. Byers and M. E. Flatte, Phys. Rev. Lett. 74, 306 

(1995); T. Martin, Phys. Lett. A 220, 137 (1996); M. P. 

Anantram and S. Datta, Phys. Rev. B 53, 16390 (1996); 

G. Deutscher and D. Feinberg, Appl. Phys. Lett. 76, 487 

(2000) ; R Melin, J. Phys.: Condens. Matter 13, 6445 

(2001) ; 255 (2001). 

[11] G. B. Lesovik, T. Martin, and G. Blatter, Eur. Phys. J. 
B 24, 287 (2001); P. Recher, E. V. Sukhorukov, and D. 
Loss, Phys. Rev. B 63, 165314 (2001); P. Samuelsson, E. 
V. Sukhorukov, and M. Biittiker, Phys. Rev. Lett. 91, 
157002 (2003); 

[12] D. Beckmann, H. B. Weber, and H. v. Lohneysen, Phys. 
Rev. Lett. 93, 197003 (2004); S. Russo, M. Kroug, T. 
M. Klapwijk, and A. F. Morpurgo, Phys. Rev. Lett. 95, 
027002 (2005); P. C. Zimansky and V. Chandrasekhar, 
Phys. Rev. Lett. 97, 237003 (2006); A. Das, Y. Ronen, M. 
Heiblum, D. Mahalu, A. V. Kretinin, and H. Shtrikman, 



Nature Comm. 3, 1165 (2012). 
[13] S. Duhot, F. Lefloch, and M. Houzet, Phys. Rev. Lett. 

102, 086804 (2009). 
[14] M. Houzet and P. Samuelsson, Phys. Rev. B82, 060517 

(2010). 

[15] N. M. Chtchelkatchev, T. I. Baturina, A. Glatz, and V. 

M. Vinokur, Phys. Rev. B 82, 024526 (2010). 
[16] B. Kaviraj, O. Coupiac, H. Courtois, and F. Lefloch, 

Phys. Rev. Lett. 107, 077005 (2011). 
[17] J. C. Cuevas and H. Pothier, Phys. Rev. B75, 174513 

(2007). 

[18] J. P. Cleuziou, W. Wernsdorfer, V. Bouchiat, T. On- 
darguhu, and M. Monthioux, Nature Nanotechnology 1, 
53 (2006). 

[19] L. G. Herrmann, F. Portier, P. Roche, A. Levy Yeyati, 

T. Kontos, and C. Strunk, Phys. Rev. Lett. 104, 026801 

(2010); P. Burset, W.J. Herrera, and A. Levy Yeyati, 

Phys. Rev. B 84, 115448 (2011) 
[20] L. Hofstetter, S. Csonka, J. Nygard, and C. 

Schonenberger, Nature (London) 461, 960 (2009); L. 

Hofstetter, S. Csonka, A. Baumgartner, G. Fiilop, S. d 

HoUosy, J. Nygard, and C. Schonenberger , Phys. Rev. 

Lett. 107, 136801 (2011); A. Das, Y. Ronen, and M. 

Heiblum, arXiv cond-mat:1205-2455. 
[21] A. Freyn, B. Dougot, D. Feinberg, and R. Melin, Phys. 

Rev. Lett. 106, 257005 (2011). 
[22] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 

(1992). 

[23] T. Jonckheere, A. Zazunov, K. V. Bayandin, V. 
Shumeiko, and T. Martin Phys. Rev. B 80, 184510 
(2009). 

[24] D. ChevaUier, J. Rech, T. Jonckheere, and T. Martin, 

Phys. Rev. B 83 125421 (2011). 
[25] L. Glazman and K. Matveev, JETP Lett. 49, 659 (1989). 
[26] J. Rech, D. ChevaUier, T. Jonckheere, and T. Martin, 

Phys. Rev. B 85, 035419 (2012). 
[27] A. Cottet, Phys. Rev. B 86, 075107 (2012); A. Cottet, 

T. Kontos, and A. Levy-Yeyati, Phys. Rev. Lett. 108, 

166803 (2012). 



